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^^ . The Hg ion plays a key role in the chemistry of dense interstellar gas clouds where star; 

,.£^ ' and planets are forming. The low temperatures and high extinctions of such clouds make 

Mh. direct observations of H;^ impossible, but lead to large abundances of H2D+ and D2H^ 

Q ' which are very useful probes of the early stages of star and planet formation. Maps o 

H . H2D+ and D2H+ pure rotational line emission toward star-forming regions show thai 

^ ' the strong deuteration of H^ is the result of near-complete molecular depletion of CNO 

bearing molecules onto grain surfaces, which quickly disappears as cores warm up aftei 

stars have formed. 

In the warmer parts of interstellar gas clouds, HJ transfers its proton to other neutral; 
such as CO and N2, leading to a rich ionic chemistry. The abundances of such species are 
useful tracers of physical conditions such as the radiation field and the electron fraction 
QQ ' Recent observations of HF line emission toward the Orion Bar imply a high electroi 

\^ . fraction, and we suggest that observations of 0H+ and H2O+ emission may be used tc 

r^ ' probe the electron density in the nuclei of external galaxies. 
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1. Introduction 

The Hj^ ion plays a central role in the physic s and chemis try of interstella: 
gas clouds where stars and planets are forming ( Watsonlll973f ). In cold parts o: 



such clouds, where temperatures are ^20K, H3 is converted into its deuteratec 
forms, which become the dominant forms of Hj^ at T^IOK. Here, H^ can bf 
studied directly by rotational spectroscopy of H2D"^ and D2H"'" in the far-infrarec 
and submillimeter. In warm interstellar clouds (T~30-100K), a direct study o: 
Hj^ is possible through vibrational spectroscopy in the mid-infrared, and thit 
technique is widely used for lines of sight through the diffuse interstellar mediun 
(see contributions by McCall, Indriolo and Geballe to this volume). To peer intc 
highly opaque star-forming regions and to map their spatial structure, emissior 
from indirect probes such as HCO^ and N2H"'" is used, which form by protor 
transfer of Hg" to CO and N2, respectively. 

The goal of this paper is to review recent work on the astrophysical use of H3 
and its deuterated forms. Since the role of Hjj' in the warm and cold regimes if 
so different, this review discusses them separately: § |2] covers the study of H2D^ 
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and D2H"^ in cold clouds, while §§ |3] and |4] treat probes of Hjj" in warm regions. 
Space limitations prevent a complete coverage of all work since the previous H^ 
conference in 2006; therefore the content of this review is necessarily biased toward 
the interests of its author. In particular, this paper does not cover the recent results 
on the ortho/para ratios of H2 and H2D"'", which may constrain the ages of pre- 
stellar cores; see the contribution by Pagani to this volume. This work is also of 
interest as it marks the first use of the deta iled state-tq-state rate coefficients for 
de-excitation of H2D"'" in collisions with H2 ( Hugo et al.ll2009l ). Collisions between 
Hg" and H2 may lead to reaction (fully elastic case), excitation (fully inelastic 
case), or both. Thus, the rates by Hugo et al. are a major step forward for the 
interpretation of interstellar H2D"'" lines, compared to the s caled radiative rates 
used previously (see discussion by IVan der Tak et all |2005[ ) . The next step up 
from Hugo's statistical approach would be full quantum scattering calculations, 
which would be a large effort because of the reactive nature of deuterated Hg" . 

Another important topic that this paper does not discuss in detail are 
observations of H2D^ in protoplanetary disks. The conditions in the midplanes of 
such disks are somewhat similar to those in pre-stellar cores: the densities are high 
enough {Z 10^ cm~'^) that ultraviolet radiation does not penetrate, leading to a low 
temperature (<10K) and a low ionization fraction. The theory of chemistry in 
the limit of complete depletion of CNO-bearing species ( Walmslev et al .112 004 ) 
was originally developed to describe pre-stellar cores, but may be even more 
applicable to protoplanetary disk midplanes. Here, all but the lightest molecular 
species freeze out onto grain surfaces, and H2D"^ is the only available ion to 
probe the ionization rate of the gas, which controls the efficiency of interaction 
with the magnetic field. Since the magneto-rotational instability is thought to 
play a major role in the dynamics of such disks, several team s have searched 
for H ^D"*" line emission. After an initial clain ied detection by ICeccarelli et al 



( 2004} ). only upper limits have been reported bv iGuilloteau et al] ((2006), Qi et al. 



(|2008[ ). IChapillon et~aD (l201lf ) and lOberg et all (|201lh . See these papers for the 
quantitative implications of the limits, and the perspectives with future telescopes. 



2. Cold gas: Spatial distribution of H2D"*" and D2H+ 



As reviewed previously (jVan der Takl 120061 ) , the astrophysical use of deuterated 



H^ is based on the fact that the formation of H2D"'" and D2H"'" from Hj^ + HD 
is strongly enhanced at low temperatures (<20K). As a result, the abundance 
of deuterated Hj^ peaks in the dense, cold concentrations of gas within clouds 
known as cloud cores. The temperature dependence is enhanced by a second 
effect: the main destroyer of deuterated H^j" is CO, which freezes out onto grain 
surfaces under the same dense and cold conditions. These two effect make H^D"*" 
and P 2H"'" ideal tracers of the initial conditions of star formation; see lSipila et al.l 
( 2010f ) for recent models of this situation. At higher temperatures (~20-30K), 
reactions wi th CH2D"'" and CoHD"*" rather than H2D^ act to enr ich molecular 
D/H ratios (jParise et al.ll2009l : ISakai et al.ll2009l : lOberg et al.ll2012h . 

The importanc e of deuterated H|^ has been realized for a long time e.g., 
Pagani et al.lll992l but whether this potential coul d be used remained unclear until 
the first detection spectra of interstellar H2D"'" ( Stark et al.lll999l : ICaselli et al 
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20031 ). Motivated by these results, ICaselli et al.l ( 20081 ) present the first survej 



of H2D"'" line emission, toward a sample of 10 starless cores and 6 protostellai 
cores. Using the Caltech Submillimeter Observatory, the H2D"'' line at 372 GHz it 
detected in 7 starless cores and in 4 protostellar cores. The column density of o- 
H2D+, derived from the 372 GHz line intensity, is also higher in starless cores thar 
in protostellar cores. The brightest H2D+ lines are detected toward the densesi 
and most centrally concentrated starless cores, where the CO depletion factoi 
(measured by the ratio of CO line emission to dust continuum emission) and th( 
deuterium fractionation (measured as the N2D"'"/N2H"'" line ratio) are also largest 

The results of this survey are consistent with the expectation that H2D"'" trace; 
very dense and cold gas, partly because the reaction of H|^ + HD proceeds mainlj 
in the forward direction, and partly because the main destroyer of H2D"*",C0 
freezes out onto grain surfaces. In particular, at gas temperatures above 15 K, low 
CO depletion factors and large abundance of negatively charged small dust grains 
or PAHs drastically reduce the deuterium fractionation to values inconsistent witl 
those observed toward pre-stellar and protostellar cores. 

Comparing their results to a chemical model, ICaselli et al.l ( 20081 ) find thai 
several parameters, including the cosmic-ray ionization rate, the depletion o: 
CO (and generally of neutral species), the volume density, the dust grain siz( 
distribution and the fraction of PAHs all affect t he H2D"'" abundance . In addition 



the ortho-to-para ratio of H2 is a key parameter ( Pagani et al.ll2009l ). The curreni 
data do not allow to disentangle these effects: observations of the P-H2D+ grounc 
state line would help to break this degeneracy, as would detailed measurements 
of the temperature s and densitie s of th e cores as a function of radius. 

The results by ICaselli et al.l ( 20081) are corroborated by mapping of H2D^ 



line emission in the Oph B2 core ( Friesen et al.l I2OIOI ). The column densitj 



distributions of both H2D"'" and N2D"'" in this core show no correlation with th( 
total H2 column density. Instead, the deuterium fractionation in the Oph B^ 
core, again measured by the N2D"'"/N2H"'" ratio, systematically decreases witl 
proximity to the embedded protostars, out to distances ^0.04pc. 

Us ing the multi-pixel CHAMP+ receiver on the APEX telescope. iParise et al 
( 201ll ) have successfully detected D2H"'" emission at 692 GHz from the H-MM] 



prestellar core located in the LDN 1688 cloud. The emission is detected or 
several pixels of the CHAMP+ array, hence extended on a scale of at least 40" 
corresponding to ~4800 AU. T his is the fir st secure detection of interstellar D2H^ 
after the tentative detection bv lVastel et a l. (2004). The column densities of H2D"^ 
and D2H+ toward this core are similar, which implies strong depletion of CO 
model calculations indicate that 90-99% of CO is frozen on dust grains. 

Other maps of H2D"'" emission will be published soon. The region NGC 1331; 
IRAS4, where H?D" ^ was first detecte d, is part of the Spectral Legacy Survej 
(SLS) at the JCMT (|Plume et al.ll2007l ). which initially covered the 330-360 GH2 



range but was later extended to 373 GHz. This survey provides images of moleculai 
line emission over a 2x2' field at 15" resolution toward four different star-forming 
environments. In the case of NGC 1333 IRAS4, the H2D"'" emission is found tc 
have a spatial distribution unlike that of other ions such as DCO'*' or N2H"'" or th( 
dust continuum, which peaks near the embedded protostars. The fact that H2D^ 
emission peaks away from heat sources confirms the idea that H2D"'" is rapidlj 
destroyed as cores warm up (E. Koumpia et al, in preparation). 
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The first mapping survey of H2D+ toward pre-stellar cores is also underway. 
Using tlie JCMT, Di Francesco et al (in preparation) liave made maps of H2D+ 
line emis sion toward a s ample of 6 pre-stellar cores, selected from the N2D"^/N2H+ 
study bv lCrapsi et al.l (J2005l ) to have high N2D"'"/N2H"'" ratios. For each of these 
cores, the H2D"'" emission is found to peak at the core center. The emission is 
extended over ~30", as seen in the maps and confirmed by the similarity of H2D~'" 
antenna temperatures between the CSO and JCMT telescopes. 



3. Warm gas: Reactive ions in star-forming regions 

Whereas large abundances of interstellar H2D^ and D2H^ require low 
temperatures ( J; 10 K), proton transfer from Hg" to other molecules also takes place 
in warm gas. In particular, reactions of interstellar Hj^ with species such as CO and 
N2 lead to HCO"*" and N2H"'" which are widely observed in the interstellar medium. 
Stable ionic species a re useful as tracers of the interaction of inters tellar gas with 
magnetic fields (e.g.. lHoude et al]l2004l : ISchmid-Burgk et al.ll2004l ). whereas ions 



which react rapidly with H2 trace other physical conditions such as the gas density 
and the local radiation field. The second part of this review paper discusses a few 
recent results in our understanding of the physics and chemistry of reactive ions 
in interstellar space. 

(a) The HnO^ puzzle 

The HIFI instrument is a heterodyne spectrometer onboard ESA's Herschel 
space observatory, which provides access to the spectral ranges between 480- 
1250 and 1410-1910 GHz at a resolutior i of R^O.l MHz. One HIFI highlight is th e 
discovery of interstellar 0H+ and H2O+ (|Gerin et al.ll2010l : IOssenkoDf et al.ll2010f ). 



which are intermediate steps in the ionic formation route of interstellar H2O. This 
route begins with charge transfer of H"*" or Hg to O, proceeds with reactions 
of H2 with O"*", OH"*" and H2O"'", and ends with the dissociative recombination 
of HsO"*" which produces H2O. Although this route dominates the formation of 
H2O in the gas phase at temperatures ^250 K, which is the bulk part of the 
neutral interstellar medium, the abundances of OH"*" and H2O"'" were predicted 
to be low because these species react rapidly with H2, which is the bulk species 
of dense interstellar gas clouds. In other words, the rate-limiting steps of ionic 
H2O formation were thought to be the formation of O"^ and the dissociative 
recombination of HaO"*". 

Observations with HIFI have revealed strong abs orption by interstellar 0H+ 
and H2O"'" on many lines of sight in our Ga laxy (Neufeld et al.l l2010al ). and 



toward nearby star burst galaxies such as M82 ( Wei£ et al.ll2010l ). NGC 253 and 
NGC 4945. The large column densities of these reactive species imply that the 
hydrogen in these clouds is neither purely in atomic form, because no OH"*" and 
H20~'' would be produced, nor in purely molecular form, because all OH"*" and 
H2O"'" would react into HsO"*" and H2O. Previous models of the structure of 
molecular clouds, where clouds are either diffuse and hydrogen atomic, or dense 
and hydrogen molecular, a re inconsistent with these observations and must be 



rejected (e.g.. lSnow fc McC all 2006). While PDR models are able to explain the 
observed amounts of OH"*" and H20~'' with diffuse gas where most hydrogen is 
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Figure 1. Spectrum of the HF J=l-0 line, observed with Herschel-HIFI toward the Orion Bar 

From: Van der Tak et al (2012). 



atomic (Gerin, this volume), the distribution of diffuse and dense gas must be 
reconsi dered, as aheady ind icated by mm-wave absorption hne studies of diffuse 
clouds (iLucas fc Lisztlll996[ ). 

A major puzzle in the study of interstellar 0H+ and H2O"'' (which I collectivelj 
call H^O"*") is the Herschel-SPIRE spectrum of the active galactic nucle us Mrk 231 
where lines of H„0+ appear in emission ( Van der Werf et al.l |2010| ). The largf 
dipole moments and the small reduced masses of HnO"*" imply high line frequencies 
and large radiative decay rates, so that excitation of their rotational levels requires 
extremely high densities and line emission is not expected to be observable. Since 
my previo us discussion of thi s topic at the SMILES conference, also organized bj 
Prof. Oka ( Van der Takll2010l ). a seco nd report of HnO"*" eni ission has appeared fo] 
the ultraluminous merger Arp 220 ( Rangwala et al.l l20 111 ). In this galaxy, botl 
the OH"*" and the H20^ lines show a P Cygni-type profile, which consists o: 
blueshifted absorption and redshifted emission, and which is characteristic of £ 
wind or an outflow. 

One hint at the solution of this "H„0"^ puzzle" may be the HF molecule. Like 
H^O"'', HF appears in a bsorption on most lines of sight through the interstellai 
medium of our Galaxy ( Sonnentrucker et al.ll2010l ). which is expected from thf 
large dipole moment and small reduced mass of t he molecule. The SPIRE 
spectrum of Mrk 231 however shows HF in emission ( Van der Werf et al.ll201[ll 
and the SPIRE spect rum of Arp 220 shows a P Cygni profile for the HF linf 
(JRangwala et al.ll201ll ). This behaviour is the same as for H^O"*", which suggests 
(but does not imply) a connection between these species. 



(6) Detection of HF emission 

In dense interstellar gas clouds, the bulk of the hydrogen is in the form of H2 
which is essentially unobservable at low temperatures. The common tracer for H; 
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clouds is the CO molecule, which is abundant and chemically stabl e, but which 
photo dissociates at H2 column densities below ~1.5xl0^^ cm~^ (|Visser et al. 



20091). While CH m ay be used for the intermediate (translucent cloud) regime 
([Sheffer et al.l 120081 ) . a proxy for H2 at very low column densities was missing 
until the Herschel mission. Since the reaction of F with H2 leading to HF is 
exothermic and HF photodissociation is slow, HF is expected to be the main 
carrier of fluorine in the gas ph ase, as long as ^ 10~^ of hydrogen is in the form 
of H2 (JNeufeld fc Wolfirell2009f ). 



The first confirmation of these predictions came from observations with 
the LWS instrument onboard the ISO satellite, which me asured absorption in 
the J=2-l line of HF toward Sgr B2 (iNeufeld et al.lll997D . The Herschel-HIFI 
instrument gives us first access to the rotational ground state of HF, and 
observations of widespread absorption in this line indicate an HF abundance of 
(l-2)xl0~^ in diffuse clouds, close to th e interstellar (and the solar) abundance 
of fluorine ( Sonnentrucker et al.l l2010f ). Absorption is al so observed towar d 
dense clouds, but the implied abundance is ~100x lower (jPhillips et al.ll2010f ). 
suggesting that depletion on grain surfaces plays a role. 

The large dipole moment of HF and the high frequency of its J=l-0 line 
imply that radiative decay to the ground state is fast, which explains why the line 
usually ap pears in absorp tion on lines o f sight through the interstellar medium 
( Neufeld et al. 2010b ; Monje et al.ll2011al ). Only very dense enviro nments such as 



the in ner envelopes of late-type stars show the HF line in emission (lAgundez et al. 



201 ll ). The only detection of HF emission from the Galactic interstellar medium 
so far is toward the Orion Bar, and it is a surpris e because the H2 densit y 
in this region is not high enough to excite the line ( Van der Tak et al.l l2012h . 
Instead, collisions with electrons appear to dominate the excitation of HF; non- 
thermal excitation mechanisms seem less likely, as detailed in § |4]of this paper. 
Line emission of HF thus appears to trace regions of molecular gas with a high 
electron density, which in the case of the Orion Bar is caused by strong ultraviolet 
irradiation by the Trapezium stars. 

Excitation by electrons may also apply to active galactic nuclei where 
the Herschel-SPIRE spe ctra show HF in emission such as Mrk 231 



( Van der Werf et al.l I2OIOI ) or as a P Cygni profile such as Arp 220 
Rangwala et al.ll201lr). In coritrast, t he Clover leaf quasar at z = 2.56 shows HF in 



pure absorption (JMonie et al.ll2011bl ). which may imply a lower electron density. 
Given the importance of electron excitation for HF, the appearance of the line in 
emission or absorption may be a measure of the local ionization rate by ultraviolet 
photons (in a starburst) or by cosmic rays (in an AGN). 

The similarity of the line profiles of HF and H„0^ in the Herschel-SPIRE 
spectra of Mrk 231 and Arp 220 may not be a coincidence. Electron impact 
excitation may also play a role for H^O"*" in these sources. Cross sections for 
the inelastic collision of the e-H^O"*" system are urgently needed to test this 
hypothesis. Weak OH"*" emission is also seen toward the Orion Bar, providing 
further evidence for this trend (Nagy et al, in prep.). 
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4. Alternative mechanisms to produce HF emission 

(a) Infrared pumping of HF 

Radiative excitation of tlie HF line may occur tlirougli tlie f=l-0 vibratior 
mode at 2.55 //m. Tliis wavelengtli is too sliort for resonances witli PAH features 
wliicli leaves dust continuum and H2 line emission as options. Infrared pumping 
is not known to play a role for other molecular species in the Orion Bar, but mos1 
of these species occur deep in the molecular cloud where near-infrared radiatior 
does not penetrate due to extinction by dust. In contrast, HF readily forms a1 
low extinction where both dust continuum and H2 line emission are strong. 

We have c a lculat ed the positions of the t)=l"0 lines of HF from the constants 
of iRam et aD (jl996l ) up to the J— 10 level. Three lines are close to a line o: 



H2: the P(l) line at 3920.312 cm'^ the R(0) line at 4000.989 cm-\ and th( 
R(l) line at 4038.962 cm^^. Their respective partners are the H2 v=2~l Q(l' 
line at 3920.053 cm"i, the H2 v=l-0 Q(7) line at 4000.075 cm-\ and the hJ 
?;=l-0 Q(6) line at 4039.507 cm~^. Especially the first match is close, which it 
encouraging because the P(l) line leads directly to the J=l level observed witl 
HIFI, and the m atching H2 line is one of the strongest H2 lines in Model 14 o: 
iBlack &: van Dishoeck (1987i ). with an intensity of ~60% of the well-known t>=l-C 
S(l) line at 2.12 //m. However, the separations of these line pairs of 20, 69 anc 
40 kms~^ exce ed the widths of the H2 pure rotational lines in the Orion Bar o: 



4-6 km s ^ (Allers et al.ll2005l ). The spatial distrib ution of the H2 pure rotati ona 



lines is similar to that of the rovibration al lines ( Van der Werf et al.l Il996l ) , sc 
that their widths are probably similar too; iTielens et al.l ( 19931 ) already rules ou1 



a large contribution by shocks to the H2 rovibrational line emission from the Orior 
Bar. 

Alternatively, the infrared pumping of HF may take place by continuun 
radiation from dust. Following ICarroU fc GoldsmithI (J1981 ). efficient pumping 
requires that ^vibe//(e ' '^ — 1) exceeds ^rotj where ^vib and A^ot are the 
vibrational and rotational Einstein A coefficients, e is the dust emissivity (assumec 
to be unity in the mid-infrared), T^ is the dust temperature, and / is the geometric 
filling factor of the dust. Since the vibrational frequency of 3961.4 cm~^ is higl 
[hv/k » Td), the expression simplifies to /e~ ' '' > Amt/A y ^h^ wh ere the righ1 
hand side is about 1/100 using Ayi]^—203s~^ from iRam et al.l ( 19961 ). 



If the dust radiation field fills the entire sky as seen from the HF cloud (/=1) 
the lower limit for the dust temperature is 1170 K. While such warm dust coulc 
exist locally in the Orion Bar, this temperature seems unreasonably high as ar 
average over a 7200 AU region. Even transiently heated mini-grains will rarelj 
get so warm. Smaller filling factors exacerbate this problem: for example 50% skj 
coverage means a dust temperature of at least 1380 K, and 10% coverage impliet 
T[i >2340K. Hence continuum radiation cannot do the pumping either. 

Finally, we consider the possibility that the excitation of HF in the Orior 
Bar is influenced by radiation from the Trapezium stars. T he brightest of thes e 
stars is 6^C Ori, which is 127" away from the HI FI position (Van LeeuwenI 120071 ) 



Assuming that both objects lie 414 pc from us (JMenten et al.ll2007[ ). the stella: 



continuum flux at the Orion Bar is 2.6 million times stronger than at the Earth 



For t he 2.5 /im flux density, we adopt the K band magnitude of 4.14 (iMuench et al 



l2002h which corresponds to 11.5 Jy at Earth and 30.3 MJy at the Bar. Dividin, 
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this flux density by Att gives an angle-averaged radiation field Jj, of 2.4x10^^^ 
ergs~^ cni~^ Hz~^ sr~^. Setting this equal to the Planck function at a radiation 
temperature Tr, we obtain rR = 209K at A = 2.5 //m, which is only slightly 
larger than the intensity due to average Galactic starlight at the Sun's location. 
The corresponding pumping rate in the HF vibrational fundamental is of order 
Avih / {^^p{hiy / kTji) — 1} = 1.8x10"^^ s~^, which is much less than the collisional 
excitation rates by H2 and electrons. We conclude that infrared pumping does not 
play a role for HF in the Orion Bar. 

(6) Formation pumping of HF 

Besides infrared pumping, an alternative non-thermal excitation mechanism 
of HF would be chemical pumping. In this scheme, the excitation of HF is out 
of thermal equilibrium because the formation and destruction rates are as fast 
as radiative decay. This mechanism is plausible for the Orion Bar because of the 
high ultraviolet radiation field, which leads to a large amount of material in the 
regime where H2 is abundant but the HF destruction rate by photodissociation 
and reaction with C"*" is of the order of the maximum formation rate. 

Assuming optically thin emission, the observed HF line flux of 7.8Kkms~^ = 
1.5xl0~^ ergcm"^ s~^ sr~^ requires a column density of Nj=i = 9.5x10^^ cm~^ 
in the J=l state of HF. To obtain this by chemical pumping, the HF formation 
rate R = kn(R2)n{F)/n{RF) must fulfil R/Aw x X^f x Nh> 9.5 x lO^^cm'^, 
where ^10 is the Einstein A coefficient of the HF J=l-0 line. The formation of 
HF is dominated b y the reaction of Ho with F, the rate coefficient of which is 



given in Table 2 of lNeufeld fc Wolfird (|2009l ) and equals k = 4.6x10-^^ cm^s-i 



at T = 80 K. In strong radiation fields, destruction of HF is dominated by 
photodissociation, at a rate of Cd = 1.2 xlO"^*^ s~^ in the standard Draine field, so 
that (d — 6xl0~^s~^ for the Orion Bar (x~50,000). The HF abundance follows 
from balancing its formation and destruction rates {(d = R): 

n{RF)Cd = kn{R2)n{F) (4.1) 

where the sum of the abundances X(HF) -|- ^(F) is constrained by the total 
amount of fluorine, 1.8 x 10~^ relative to H. The product R/Aiq x X(HF) reaches 
its maximum of 4.5x10"^^ for very high densities, leading to a maximum 
Nj=i = 4.5 X 10^^ cm~^ when assuming a total gas column density of 10^^ cm~^. 
For gas densities of 10^ cm~^ this drops by a factor of two, for a gas density 
of 10^ cm~^ by another factor of ten. Higher temperatures lower the required 
densities approximately proportionally. 

The computation above assumes that every molecule formation leads to an 
excitation of the J=l state. This premise is reasonable, if not intuitive. The 
formation of HF releases 1.36 eV in each reaction, which is well above the energy 
of the u=l state, so that all rotational levels in the vibrational ground state will 
be about equally populated. The relative fraction of molecules formed in J=0 will 
be very small, and the radiative decay of all other molecules passes through the 
J = 1 state so that they contribute to the observed J=l-0 emission. 

We have calculated the formation rate of H F in a self-consis tent model 
of fluorine chemistry using the KOSMA-r code ( RoUig et al.l l2006l ) . The rate 



equations are solved for a spherical clump of gas and dust with a density of 
nif=10^cm~^ at the outer radius (^^=0) and increasing to 1.1x10^ cm~^ at the 
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Figure 2. Chemical model of the Orion Bar as a function of Ay , using the UMIST reactioi 
rates (dashed lines) or the H2 + F rate from iNeufeld et alj (|2005l ) (solid lines). Top: formatioi 
rate of HF due to reactions of F with H2 in its w=0 and v=\ states. Middle: absolute densitiei 
of H2 and HF The apparent dip in the HF abundance near Ay = 3 when using UMIST ratei 
is due to incorrect temperature dependence of the H2 + F reaction rate. Bottom: calculatec 
temperatures of gas and dust. The distribution of HF is unlike that of OH"*" and H20^, whicl 
for this X peak at 0.05 < Ay < 1 and disappear at Ay > 1 where H2 and HF are abundant (M 
RoUig, priv. comm.). The observational similarity between H„0^ and HF is thus likely related t( 
the non-uniform structure of the ISM and similar excitation requirements, rather than chemica 
similarities. 
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clu mp center, so that the t o tal cl ump mass is IMq. The dust model is number 
7 of IWeingartner &: Draine ( 200lh. eau i valeii t to -Ry=3.1, and the formation of 
H2 is according to ICazaux fc Tielensl (l2002h with formation on PAH surfaces 
suppressed. The total abundance of fluor ine in the Orion region is taken from 
6.68 xlO~^ (ISim on-Diaz fc Stasihska"201lV which is lower than the average value 
from lNeufeldfc Wolfire (2009 ) by about a factor of 3. 

Figure [2] shows the results of our calculations, which consider both ground- 
state and vibrationally excited H2. Collision of F with H2 (w>0) dominates the 
formation of HF until Ay=0.0036, after which point ground-state H2 takes over. 
The full model confirms the numbers from the estimate above with a formation 
rate R of about 10~^ s~^, and approximately half of the fluorine in the form of HF 
in the reactive zone between Ay ~ 0.01 ... 2 containing most of the mass. Hence, 
we find that chemical pumping of the HF 1—0 transition is effective in the Orion 
Bar, but that estimates based on current molecular data fall short by a factor of a 
few. This could be due to uncertainties in the total fluorine abundance, or higher 
temperatures and radiation fields in the interclump gas. 



5. Conclusions and future work 



The study of interstellar deuterated H-|" has come a long way in the past 16 years, 
and the conferences that Prof. Oka has organized every six years provide clear 
milestones of this progress. At the 88* birthday of the discovery of H j" in 2000, 
the first detection of interstellar H2D"'" was announced ( Stark et al.lll999l), which 
followed shortly after the discovery of interstellar Hj^ itself ( Geballe &: Okalll996f ). 
The 94* bir thday saw the firs t report of strong H2D"'" emission fr om cold pre- 
stellar cores (ICaselh et al.ll2003f ). the first H2D+ map of LDN 1 544 bvlVastel et al. 
(|2006f ) and the first tentative detection of interstellar D2H+ (jVastel et al.ll2004^ . 
This y ear, at the 100* birthday, we are seeing the confirm ation of interstella r 
D2H+ (iParise et al.ll201ll ). the first survey of H2D+ emission (ICaselh et al.ll2008f ). 
and mapping surveys of H2D"'" are ongoing. 

The natural question to ask now is: What will the 106*^^ birthday in 6 years' 
time bring? The most likely development in the next years are high-resolution 
maps of pre-stellar H2D"'' with the ALMA interferometer, which is now about 
half-way completed and has started its Early Science program. In fact, the tuning 
range of the ALMA receivers was especially extended to include the 372 GHz 
line of 0-H2D"'", following the CSO detection of strong emission in this line in the 
pre-stellar core LDN 1544. The ALMA telescope will also likely make the first 
detection of H2D"'" in a protoplanetary disk, and hence constrain the ionization 
fraction of the disk midplane. Another possibility is a survey of the P-D2H"'' line 
at 692 GHz in pre-stellar cores with ALMA, although the emission is probably 
quite weak in many cases. Third is the new capability of the GREAT receiver on 
the SOFIA airborne observatory, which gives access to the P-H2D+ ground state 
line near 1370 GHz. The science operations of SOFIA have started in 2011, and 
GREAT is among the first-light instruments. 

The GREAT instrument may also be used for searches for the ground-state 
line of o-D2H^" at 1477 GHz. A search for thi s line with HIFI in the framework 
of the CHESS program ( Ceccarelli et al.ll2010l ) was unsuccessful (C. Vastel, priv. 
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comm.). Both the P-H2D+ 1370 GHz hne and the 0-D2H+ 1477 GHz line an 
unlikely to appear in emission: their upper level energies of ~70 K and radiative 
decay rates of ~3xl0~^s~^ imply that their excitation requires higher densities 
and temperatures than occur in regions where the species are abundant. However 
it is possible that the GREAT observations will reveal absorption by p-H2D"^ 
and 0-D2H"'" near 1400 GHz. Appearance in absorption requires a source o: 
background continuum, which the warm central regions of young protostellai 
cores may provide, as well as sufficient abundance of H2D'*' amd D2H"'" in th( 
cold outer layers of these cores. Such conditions may well occur along some lines 
of sigh t, as the example of H2D"^ absorption at 372 GHz toward IRAS 16291; 
shows ( Stark et al.ll2004l ). These observations would be a new astronomical use o: 



deuterated H3 to probe physical conditions in star-forming regions. 

In warm regions of the interstellar medium, hydride molecules such as HI 
and reactive ions such as OH"*" and H20~'" will continue to play an importani 
role in characterizing the physical conditions of the gas. Mapping the spatia 
distribution of HF in the Orion Bar and other photon-dominated regions will 
Herschel-HIFI will help to clarify the contributions by atomic and molecular layers 
to the HF emission. The Herschel-SPIRE spectra of other galactic nuclei will help 
us to understand the origin of OH"*" and H20''~ emission by probing a range o: 
environments. After the end of the Herschel mis sion in R:i March 2013, ground- 
based observations of OH"*" such as pioneered bv IWvrowski et aP ( 2010f ) will b( 



invaluable as probes of interstellar physics and chemistry, suppo rted by othe] 
react ive ions that are observable from the ground s uch as SH"^ (iMenten et al 



201 Ih . C0+ (IStauber fc Brudererl I2OO9I ) and HOC+ fiRizzo et al.ll2003f ). Clearly 



the use of deuterated H^ and other molecular ions to probe star-forming regions 
has a bright future! 
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